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Abstract 

We study the dynamical evolution of a phase-transition-induced collapse neutron star to a hybrid 
star, which consists of a mixture of hadronic matter and strange quark matter. The collapse 
is triggered by a sudden change of equation of state, which result in a large amplitude stellar 
oscillation. The evolution of the system is simulated by using a 3D Newtonian hydrodynamic 
code with a high resolution shock capture scheme. We find that both the temperature and the 
density at the neutrinosphere are oscillating with acoustic frequency. However, they are nearly 
180° out of phase. Consequently, extremely intense, pulsating neutrino/antineutrino fluxes will 
be emitted periodically. Since the energy and density of neutrinos at the peaks of the pulsating 
fluxes are much higher than the non-oscillating case, the electron/positron pair creation rate can 
be enhanced dramatically. Some mass layers on the stellar surface can be ejected by absorbing 
energy of neutrinos and pairs. These mass ejecta can be further accelerated to relativistic speeds 
by absorbing electron/positron pairs, created by the neutrino and antineutrino annihilation outside 
the stellar surface. The possible connection between this process and the cosmological Gamma-ray 
Bursts is discussed. 

Keywords: dense matter : stars: neutron stars : stellar oscillations : phase transition- quark 
stars : gamma ray bursts. 
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I. INTRODUCTION 



Gamma-ray bursts (GRBs) are cosmic gamma ray emissions with typical fluxes of the 
order of 10~^ to 5 x 10~^ergs cm~^ with the rise time as low as 10~^ s and the duration of 
bursts from 10"^ to 10^ s. The distribution of the bursts is isotropic and they are believed 
to have a cosmological origin, with the observations showing that GRBs originate at extra- 
galactic distances. The large inferred distances imply isotropic energy losses as large as 
3 X 10^3 ergs for GRB 971214 and 3.4 x 10^^ ergs for GRB 990123. For detailed reviews 
of GRBs properties see 111 and respectively. The widely accepted interpretation of the 
phenomenology of GRBs is that the observable effects are due to the dissipation of the 
cinetic energy of a relativistically expanding fireball, whose primal cause is not yet known 



II. 

One fundamental question related to GRBs is how many intrinsically different categories 
they have. Each type may correspond to an intrinsically different type of progenitor, as 
well as to a different type of central engine. From the GRB sample collected by BATSE, 

n n 

a clear bimodal distribution of bursts was identified [2|, |3|. Two criteria have been used 
to classify the bursts. The primary criterion is duration, and a separation line of 2 s has 
been adopted to separate the double-hump duration distribution of the BATSE bursts. The 
second criterion is the hardness-usually denoted as the hardness ratio within the two energy 
bands of the detector. On average, short GRBs are harder, while long GRBs are softer. 
Hence the observations show the existence of two classes of GRBs, long-soft, and short- 
hard, respectively. The short GRBs are hard because of the harder low-energy spectral 
index of the GRB spectral function j4|. More interestingly, short GRB spectra are broadly 
similar to those of long GRBs, if only the first 2 seconds of data of the long GRBs are 
taken into account. Afterglows observation shed some light into the physical nature of these 
two types of GRBs. The host galaxies of long GRBs a.e exclusively sta. fotnnng galaxies, 
predominantly irregular dwarf galaxies [5:]. The observations suggest that most, if not all 
of the long GRBs are produced during the core collapse of massive stars, called coUapsars, 
as has been also suggested theoretically ^, Q, [s]. Both the observations and the theoretical 
models also support the idea that long GRBs are associated with supernova explosions [9|. 

The observations of SWIFT and of HETE led to a very different picture for the short 
GRBs [l, 2|. The basic result of all these observations is that short GRBs are intrinsically 
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different from the long GRBs. Most of the short GRBs are found at the outskirts of elhptical 
galaxies, where the star forming rate is very low. Even if they can be associated with star 
forming regions, they are rather far away from the active zones jl^. In several cases a 
robust host galaxy has not been identified, but the host galaxy is of early type. Deep 
supernova searches have been performed, but with negative results |j2j]- This suggests that 
short GRBs may be associated with mergers of compact objects, like neutron star-neutron 
star or neutron star-black hole mergers, white-dwarf-neutron star or black hole mergers etc. 



11 



12 



14 



15|. 



The possibility that the conversion of neutron stars to strange quark stars may be the 
energy source for the cosmological gamma ray bursts was suggested in • Neutron stars in 
low-mass X-ray binaries can accrete sufficient mass to undergo a phase transition to become 
strange stars. At the moment of its birth the strange star is very hot, with an interior 
temperature of around ~ 10^^ K. By approximating the strange matter by a free Fermi gas, 
the thermal energy of the star is Eth ~ 5 x 10^^ {p/ po)"^^^ R^T^i ergs, where p is the average 
mass density, po = 2.8 x 10^'' g/cm^ is the nuclear density, Rq is the radius of the star in units 
of 10^ cm, and Tn is the temperature in units of 10^^ K. For p = 8po, R% = 1, Tn = 1.5 the 
thermal energy of the newly formed quark star is Eth ~ 5 x lO^^ergs. In the original model 
of 13] it was assumed that the star would cool by emission of neutrinos and antineutrinos, 
and that the neutrino pair annihilation process uu —>■ e~e+ operates near the strange star 
surface. The energy deposited due to this process is Ei^2x 10^^ (To/lO^fsT) ergs ^ 10^^ 
ergs, where Tq is the initial temperature. The time scale for the deposition is around 1 s. On 
the other hand, the processes n + Ue —>■ p + e~ and p + z/g — > n + e"*" play an important role 
in the energy deposition. The integrated neutrino optical depth due to all these processes is 
T ^ 4.5 X 10~'^pi{^T^i, where pu is the mass density in units of 10^^ g/cm^. The deposition 
energy can be estimated to be E2 ~ Eth (1 — ^ 2 x 10^^ ergs. Here the value of the 
neutron drip density, pu = 4.3 has been used, and it has been assumed that all the thermal 
energy of the star is lost in neutrinos. The process 77 —>■ e~e~^ inevitably leads to the 
creation of a fireball, and this fireball will expand outward. The expanding shell interacts 
with the surrounding interstellar medium, and its kinetic energy is radiated through non- 
thermal processes in shocks. However, this model did not consider any internal dynamics, 
e.g. heat transport, viscous damping, shock dissipation etc, which can affect the neutrino 
emission dramatically. In fact the numerical simulations indicate that the temperature on 
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the neutrinosphere is rapidly changing with time, steady neutrino emission is not possible. 

The presence of oscillations of the resulting quark star produced by the phase transition 
induced collapse of a neutron star is one of the most intriguin g fe atures of the simulations 



performed with our Newtonian numerical code introduced in 



17|. Recently, the collapse 



process with a conformally flat approximation to general relativity was also simulated in 



18( 1 . The works of |17l . Il8| focus on the gravitational wave signals emitted by the collapse 
process. 

It is the purpose of the present paper to consider another important implication of this 
result, namely, the effect of the oscillations of the newly formed quark star on the neutrino 
emission. The oscillations can enhance the neutrino emission rate in a pulsating manner, 
and the neutrinos are emitted in a much shorter time scale. Therefore the neutron-quark 
phase transition in compact objects may be the energy source of GRBs. 

This paper is organized as follows. The phase transition process from neutron stars to 
hybrid stars, which consists of a mixture of strange quark matter and hadronic matter, is 
summarized in Section II. We describe our numerical code, which is used to simulate the 
dynamical evolution of star after phase transition in Section III. In Section IV, we calculate 
the neutrino and antineutrino emission from the neutrinosphere and the electron/positron 
pair creation rate due to neutrino and antineutrino annihilation process. In Section V, we 
calculate mass ejection from stellar surface by absorbing neutrinos/pairs and their subse- 
quent acceleration by the pairs outside the star. In Section VI we apply our model to GRBs. 
Finally a brief summary and discussion is presented in Section VII. 

II. DESCRIPTION OF THE PHASE TRANSITION 

The quark structure of the nucleons, suggested by quantum chromodynamics, indicates 
the possibility of a hadron-quark phase transition at high densities and/or temperatures, as 



suggested by 



19|, 



2Q|, 



2l| . Theories of the strong interaction, like, for example, the quark 
bag models, suppose that breaking of physical vacuum takes place inside hadrons. If the 
hypothesis of the quark matter is true, then some of neutron stars could actually be strange 



stars, built entirely of strange matter 22|, |23|, l3l|, l32|, l33 |. 



The central density of compact stellar objects may reach values of up to ten times nuclear- 
matter saturation density, and therefore a phase transition to deconfined quark matter, or 
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pion and kaon condensates, should take place at least in the central region of the neutron 
stars. In the case of the transition to quark matter, in addition to a phase of unpaired normal 
quark matter, present at low densities, several superconducting phases, such as the two- 
flavor color superconductor (2SC) phase or the gapless Color-flavor-lock (CFL) phase 
can also occur at the large baryon densities reached at the central regions of a compact star. 
The transition from the hadronic to the quark phase could proceed in two steps. In the first 
step, a transition from hadronic matter to normal quark matter or to a 2SC phase takes 
place, due to the increase of the baryonic density at the center of the star. This increase in 
density may be due to mass accretion from the fall-back material or rapidly spin-down of the 
star. The newly formed hybrid or quark star, containing some 2SC quark matter, can become 



meta-stable, and decay into a star containing a CFL phase [2j]. Pure hadronic compact 
stars above a threshold value of their mass are metastable. The metastability of 
hadronic stars originates from the finite size effects in the formation process of 
the first strange quark matter drop in the hadronic environment |25| . 

A phase transition between the hadronic and quark phase occurs when the pressures and 
the chemical potentials in the two phases are equal, Ph = Pq, fJ^h = f^q, where Ph, fih and Pq, 
fiq are the pressures and the chemical potentials in the hadron and quark phase, respectively. 
In the present Section, unless otherwise explicitly specified, we use the natural 
system of units with h = c = 1. If the transition pressure is less than that existing at 
the center of the compact object, the transition can occur. For finite values of the surface 
tension, complicated structures can develop in the mixed phase, such as drops, rods, and 



slabs [28|. The formation of a mixed phase is the result of two competing processes: the size 
of the barrier that the system has to overcome in order to form a structure, and the size of 
the perturbation of the system. One method to see if the phase transition can proceed or 
not is to compare the temperature reached by the system immediately after the conversion 
with the height of the barrier. If the temperature is not much lower than the height of 
the barrier, the structure formation proceeds thermally, and it is very rapid [29!,130[. If the 
temperature is low, new structures can form only via quantum nucleation, which is a very 
slow process {2^. 

The thermal nucleation rate of quark drops can be estimated in the frame- 
work of the nucleation theory as Rnuci ~ aexp {—Wc/T), where the prefactor a (the 
product of the dynamical prefactor and of the statistical prefactor) is the prod- 
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uct of a density and a growth factor, Wc = W (Re) represents the work needed to 
form the smallest drop capable of growing, and T is the temperature. Wc corre- 
sponds to the maximum of the free energy of the drop in the new phase. The 
dynamical prefactor for the nucleation rate of bubbles or droplets in first-order 
phase transitions for the case where both viscous damping and thermal dissipa- 



tion are significant has been obtained in [261]. This formalism was applied for 



the study of the nucleation of quark-gluon plasma from hadronic matter in [27| • 
By taking into account the explicit forms of the temperature dependent quark 
matter equation of state, of the dissipative factors and of the statistical prefac- 
tor, one obtains a complicated dependence of the prefactor on the temperature. 
However, the prefactor in this model is based on poorly known physical param- 
eters (like, for example, the shear viscosity of the hadronic phase). In order 
to obtain some qualitative estimates of the nucleation rate we use the thermo- 



dynamical approach developed in [291 . l30t |. The form of the prefactor can be 
obtained from general thermodynamical considerations as a ~ T^, and it is very 
little dependent on some particular choices (for example in the low temperature 
limit the prefactor in the expression for Rnud can be replaced by the baryon 
chemical potential). This form for the prefactor does not include any kinemat- 
ics, e.g. the microscopic processes required to transform a gas of hadrons into 
a gas of quarks. In the expression of the nucleation rate the dominant term is 
the exponential. The free energy is given hy W = —AuR^/S.P/'i + A'KaR^ + 87C'jR + NgA/j,, 
where AP = Pq — Ph is the pressure difference, a = (Xg + is the surface tension, 7 = 7^ — 7/1 
is the curvature energy density and A/i = fiq — fit is the difference in the chemical potential. 



24 



29|, 



3G, 



3^. 



Nq is the total baryon number in the quark drop 

The free energy has a maximum at the critical radius Rc = a (1 + VlTb) /C, where 
C = AP — UqAfi and b = 2'-fC/a'^. The corresponding free energy is 



Wr = Svra 



l + {l + bf^ + 3b/2 /3C''. (1) 



The number N of drops of the new phase formed inside the old phase in a volume in a 
time t is given by = RnudVt. Let A be the spacing between two drops in the mixed phase. 
The number of drops in a volume V is given by V/ X^. The number of drops produced while 
the front moves over a distance A must be of the order of the drops that are present in the 
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mixed phase, RnudVt = RnuciS (X/v) '>V/\^ = S/X^, where v is the velocity of the front 
and 5* its surface area. Therefore in order for thermal nucleation take place the condition 
Wc/T < ln(T^A^/f) must be satisfied. The process of absorption of a hadron into a pure 
quark matter phase can also be described phenomenologically as the fusion of a small drop 
of quarks growing into a much larger drop. 

The nature of the conversion process from the neutron to the quark matter 
is not yet fully understood, and there is no clear theoretical evidence if it is a 
deflagration or a detonation. Indeed, for realistic equations of state of quark and 
neutron matter and when matter is assumed to be in (3 equilibrium, detonation 



is difficult to achieve 



35 



36l |. However, it would be still possible to obtain a 



detonation if the matter immediately after the front is not yet in /5-equilibrium 
|24l|. It is also important to note that neutrino trapping delays j3 stability. 
According to the analysis of {2^, when the drop starts expanding, the process 
of conversion can be extremely fast, within the layer in which deconfinement is 
energetically favorable, even in the absence of the weak processes. In this case, 
the conversion front moves at the velocity of the deflagration front Vdf, which 
approaches the velocity of sound. As the conversion layer moves outward, the 
front enters the region of mixed phase, where Vdf decreases until it vanishes at 
the low density boundary of the mixed phase. 

If the conversion is indeed extremely fast, taking place at speeds close to the 
speed of sound, then the typical time scale for the transition can be estimated 
as Ttr = R/cs, where R is the radius of the neutron star and Cs is the speed of 
the sound. A simple phenomenological model for the evolution of the quark phase can be 
obtained by assuming the relation dr/dt = (r — Rc) /rtr {34 ], which gives for the transition 
time scale T^j. from a microscopic quark drop to a quark matter distribution of a macroscopic 
size the expression 



^ 10-'N;'/'Re In ^ s, (2) 



R_ 

iq 

where Rq is the neutron star radius in units of 10^ cm, Rg ~ SOORc is the initial size of the 
quark drop \m\, and Ng, which could be as large as 10"^^ |37l], is the number of quark drops 
inside the core of the neutron star. In this case R in Eq. ([2]) is replaced by R/Ny^. Therefore 
the phase transition may take place in a time scale much shorter than submillisecond. 
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III. SIMULATION OF THE PHASE-TRANSITION INDUCED COLLAPSE 



In this paper we focus on studying the phase-transition induced collapse of neutron star 
to a hybrid quark star, which consists of a mixture of strange quark matter and hadronic 
matter, and want to demonstrate that this process can produce relativistic ejecta, which 
could be a mechanism for GRBs. In order to avoid other complications we would like to 
simulate a non-rotating, phase- induced collapse of a compact object, with a mixed phase 
of quark matter and nuclear matter. In this model gravitational radiation will not be 
emitted. Hence we can focus on studying the neutrino and pair emissions and subsequently 
gamma-ray emission via interaction between the interstellar medium from this system. In 
our simulations, we will not simulate the phase transition process. Instead, we assume that 
a fast phase transition has happened (e.g., via a detonation mode) so that the initial neutron 
star has converted to a quark star in a timescale shorter than the dynamical timescale of 
the system. We assume that the normal matter inside the initial neutron star has suddenly 
changed to quark matter at t=0. This is achieved by changing the EOS at t=0 after the 
initial hydrostatic equilibrium neutron star has been constructed. We then simulate the 
resulting dynamics of the system triggered by the collapse. 



A. Description of the numerical code 

Our numerical code is based on the three-dimensional numerical simulation in Newtonian 
hydrodynamics and gravity. The quark matter of the mixed phase is described by the MIT 
bag model and the normal nuclear matter is described by an ideal fluid EOS. This code has 
been used to study the gravitational wave emission from the phase-induced collapse of the 
neutron stars • Here we briefly summarize the main equations and the numerical scheme 



involved. A detailed discussion can be found in [17 1. 

The system of equations describing the non-viscous Newtonian fluid flow is given by 

|^ + V-(pv) = 0, (3) 

-M + V.(p..v) + - = -p-, (4) 
dr 

— + V-((r + P)v) = -pvV$, (5) 
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where p is the mass density of the fluid, v is the velocity with Cartesian components Vi 
{i = 1, 2, 3), P is the fluid pressure, $ is the Newtonian potential and r is the total energy 
density, r = pe + pv^/2, and e is the internal energy per unit mass of the fluid, respectively. 
The Newtonian potential $ is obtained by solving the Poisson equation, = AirGp. The 
system is completed by specifying an equation of state P = P{p, e). 

The above hydrodynamics Eqs. ([3])- ([5]) can be rewritten in the so-called flux-conservative 
form, which can be solved numerically using quite standard high resolution shock capturing 
(HRSC) schemes. A HRSC scheme, using either exact or approximate Riemann solvers, with 
the characteristic flelds (eigenvalues) of the system, obtains the solution of a local Riemann 
problem at every cell interface of a flnite-differencing grid. Such schemes have the ability 
to resolve discontinuities in the solution (e.g., shock waves) by construction. Moreover, 
they have high accuracy in regions where the fluid flow is smooth. In general, integrating 
the hydrodynamics equations by a HRSC scheme involves the choices of an appropriate 
numerical flux formula and a reconstruction method of the state variables (p, e, v) for solving 
the Riemann problems at the cell interfaces. In our code, we use the Roe's approximate 
-liemann solver ^i8] for the numerical fluxes and the third-order piecewise parabolic method 



391] for the reconstruction. For the temporal discretization, we use a basic two-step method 



to achieve second-order accuracy in time. 

In the simulations, we introduced a very low density atmosphere outside the star. The 
"artiflcial" atmosphere is not physical, but it is important for the stability of the hydrody- 
namical code. This is due to the problem that the hydrodynamical codes cannot in general 
handle vacuum regions where the density is zero. In order to avoid a signiflcant influence 
of the atmosphere on the dynamics of the physical system, it is necessary to choose the 
density of the atmosphere patm to be much smaller than the density scale of interest. For 
the results reported in Section 3.4, we set Patmto be 3 x 10^ g/cm^. The effects of different 
atmospherical values have been compared in 4G|. 



B. Equation of state 

The equation of state (EOS) for neutron stars is highly uncertain. We could try all 
possible existing realistic EOS in our study. However, the main purpose of this paper is to 
demonstrate that during the phase-transition induced collapse of a neutron star extremely 
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intense, pulsating and very high energy neutrinos can be emitted. The effect is governed 
mainly by the amount of pressure reduction after the phase transition as compared to the 
initial neutron star model. For simplicity we will use a polytropic EOS for the initial 
equilibrium neutron star: 

P = kop'\ (6) 

where and Fq are constants. On the initial time slice, we also need to specify the specific 
internal energy e. For the polytropic EOS, the thermodynamically consistent e is given by 

^° p^^-\ (7) 



Fo-1 

Note that the pressure in Eq. ([6]) can also be written as 

P = (Fo - l)pe. (8) 

In order to describe the physical properties of the neutron stars after the 
phase transition, we consider that the star can be divided in three regions. At 
the center of the star, where p > pg, we have a pure quark core (Region I), which 
is described by the MIT bag model equation of state, so that the pressure P is 
given by 

P = P^ = l{p + pe-4B),p>p„ (9) 

where B is the bag constant, and pg is a critical density for which all the hadrons 
are deconfined into quarks. It should be noticed that Pg is not in the usual 
form of P = (ptot — 4P)/3, where ptot is the (rest frame) total energy density. It is 
because in our Newtonian simulations, we use the rest mass density p and specific 
internal energy e as fundamental variables in the hydrodynamics equations. We 
assume that the quark core is absolutely stable. The quark core is surrounded 
by a mixed phase of quark and nuclear matter (Region II) that can exist if the 
density of the region is higher than a certain critical value ptr (quark seeds can 
spontaneously produce everywhere when p > ptr)- Explicitly, the pressure in the 
mixed phase is given by 

P = aPq + (1 - a)Pn, pg> p> ptr (lO) 

where 

P, = (F„ - l)p6, (11) 
11 



and 



a 



(P - Ptr)/ {Pq - Ptr), for Ptr < P < Pq, 



(12) 



1, for < p, 

I — I 

is defined to be the scale factor of the mixed phase [17]. r„ is not necessarily 
equal to Tq. Finally, we have a normal nuclear matter region (Region III), 
extending from p < ptr to the surface of the star, so that in this region 



P = Pn, for p < ptr 



(13) 



A more detailed discussion about such hybrid quark stars can be found in pj 

The total energy density ptot, which includes the rest mass contribution, is decomposed 
as ptot = P + pe- We choose r„ < Tq in our simulations to take into account the possibility 
that the nuclear matter may not be stable during the phase transition process, and hence 
some quark seeds could appear inside the nuclear matter, or the convection, which can occur 
during the phase transition process, can mix some quark matter with the nuclear matter. In 
the presence of the quark seeds in the nuclear matter, the effective adiabatic index will be 
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43 



4^. For 



reduced. The possible values of 5^^^ range from 145 MeV to 190 MeV 41 
p > Pg, the quarks will be deconfined from nucleons. The value of Pg is model dependent; 

46|, where pnuc = 2.8 x lO^'' g cm~^ is the nuclear 



it could range from 4 to 8 pnuc. l4j, l45 
density. 

There are two issues regarding our EOS model to be addressed: (1) It should be noted 
that, for simplicity, we do not include the change in the internal energy from the phase 
transition when setting the initial data for the collapse. The binding energy released in 
the phase transition effectively leads to a slightly harder EOS due to the thermal pressure. 
This could be modeled by using a larger value of r„ (comparing to the one we used in this 
work). However, as long as the EOS after the phase transition is softer than that of the 
initial neutron star, the star will still collapse and stellar pulsations will be triggered. (2) 
Furthermore, the parameter ptr in our EOS model should be considered as the density below 
which the matter is dominated by hadronic matter. It is noted that strange quark matter 
(if it is more stable) cannot convert back to hadronic matter by decreasing the density. 
However, when a fluid element originally in the mixed phase moves to the lower density 
region (p < ptr), the fluid element will mix with a large amount of hadronic matter. Hence, 
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FIG. 1: Comparison of the equation of state of the mixed quark-hadron phase proposed in 28|] for 
K = 240 MeV and effective mass m* /m = 0.78 (dotted curve), and Eq. (jlOp . the equation of state 
if the mixed phase used in the present simulations (sohd curve). 

we approximate that the pressure in the outer region of the star is dominated by the hadronic 
part, which is modeled by an ideal gas. 

In the simulations, we choose Fq = 2, r„ = 1.85, B^/^ = 160 MeV and = 9pnuc. The 
transition density ptr is defined to be at the point where Pg vanishes initially. 

In Eq. (fTOj) . we have used a very simple linear combination of strange quark matter EOS 
and hadronic matter to represent the EOS of mixed phase, which consists of a mixture of 
strange quark drops and hadronic matter. In fact, the properties of a mixed quark-hadron 
phase and its implications for hybrid star structure were considered in . In Fig. [1] we 



compare the EOS of the mixed phase obtained in [28|] for a compression modulus K = 240 
MeV and an effective mass at saturation density of m * /m = 0.78 with Eq. (fTOj) used in 
the present simulation. We can see that these EOSs are quite close to each other. Although 
we use Eq. ( ITO|) in our simulation solely based on simplicity, we believe that the simulation 
results will not be changed qualitatively if we replace Eq. ( ITOl) by the EOS given in |28|. 

In the mixed phase we have assumed that the effective bag constant i? is a 
linearly dependent function of the density. The equation of state used by us 



reproduces quite well the EOS proposed in [28| to describe the mixed quark- 
hadron phase. It is important to note that the main purpose of our paper 
is to focus on the study of the dynamical response of the star after a sudden 
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phase transition, no matter what the fine details of the phase transition are. A 
full and exact description of the phase transition would require a very precise 
knowledge of the physical parameters describing both quark and hadronic mat- 
ter. However, taking into account the uncertainties in our present knowledge of 
the properties of matter at high densities, our investigations could certainly give 
at least a qualitative picture of the astrophysical implications of hadron-quark 
phase transitions in compact stars. 



C. The Neutrinosphere 

For a new bom compact star, the internal temperature is so high that neutrinos will be 
.apped .s,de .he .a. fo. aUea. a few seconds ,cf. H fo. a ,en.a. review). However, 
neutrinos very near the surface of the star can still escape from the star because the optical 
depth near the stellar surface is low. Quantitatively we can define a surface called the 



neutrinosphere with a radius as follows, e.g. |47l. l48| 



oo 



^eff = I drKeffir) = 1, (14) 

where the effective optical depth, Te// is defined as the inverse mean free path and the 
effective opacity, ^e// is given by 

(.e//)(r) = 1.202 xlO-Vo(r)(^)' (15) 

It is clear that this surface is a function of the temperature and of the density. 



D. Simulation Results 

The total time span of the simulations is ~5 ms, and the time step is 3.7 x 10~^ ms. For 
all the simulations we report in this paper, the grid spacing is set to be dx = 0.28 km and 
the outer boundary of the computational domain is at 27.5 km, which is about two times 
the stellar radius of our models. With the grid resolution we used for the simulations, we see 
that numerical damping becomes important after about 3 ms. We shall thus only present 
the numerical results up to 3 ms. During the simulations, a low-density atmosphere is added 
outside the neutron star. The density and temperature of the atmosphere are 3 x lO^g/cm^ 
and 0.003 MeV respectively. 
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FIG. 2: Variation of the central density with respect to the time for M = I.SSMq, (left upper 
figure), M = l.TM©, (right upper figure), M = 1.8Mq, (left lower figure), M = LOM©, (right 
lower figure). 



Fig. [2] shows the oscillations of the central density for stars with M = 1.55, 1.7, 1.8 and 
1.9Mq respectively. 

We can see that the oscillation period gradually decreases from ~ 0.32 ms to ~ 0.29 ms, 
when the mass of the star increases. It is because the radial oscillation is basically acoustic 
oscillation (cf. Lin et al. 2006) therefore frequency essentially increases with the average 

1 /2 

density. In fact the frequency is roughly scaled according to po • Hence, the oscillation 
period is smaller for more massive stars because they have higher density. We can also 
see that the oscillation amplitudes decrease slightly after 2 ms. We note that the damping 
is due partly to finite-differencing errors and physical effects, such as mass-shedding near 
the surface. One might worry that the large amplitude oscillations seen in the collapse 
simulations were numerical artifacts due to numerical errors or instability. 

In order to demonstrate that the oscillations were in fact triggered by the phase transition, 
we show in Fig. [3] the evolution of the central density (normalized by its initial value) for an 
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FIG. 3: Evolution of the central density (normalized by its initial value) for an unperturbed 
equilibrium neutron star with M = 1.8Mq. 

unperturbed equilibrium neutron star model with M = I.SMq. We see that the amplitude 
of the oscillations of the star triggered by finite-differencing errors is much smaller than that 
seen in the collapse models. Furthermore, there is no obvious periodic features in Fig. [31 

With the simulated density and temperature profiles on a given time slice, we can calculate 
the position of the neutrinosphere R^, (see Eq. (1141) ) with a trial- and-error method. Since 
Ru is a function of both temperature and density (cf. previous Section), it also oscillates 
with the same period as the central density. Figs. H] and Figs. [5] show the temperatures and 
the densities at the neutrinosphere as a function of time for neutron stars with masses of 
1.55, 1.7, 1.8 and 1.9 Mq, respectively. 

We can see that both the temperature and the density at R,, are pulsating, with the same 
period as that of the central density, but they are almost 180° out of phase to each other. 

To show the origin of the pulse like temperature and density evolution, we pick 5 time 
slices from t = to t = 0.16 ms, and we focus on the temperature and density evolution at 
the neutrinosphere. Figs. E] show the time evolution of the temperature and of the density 
profile of a star with 1.55M0 respectively. 

We find that the core temperature is rising immediately after the phase transition, which 
starts at t = 0, and the heat is moving outward. The neutrinosphere is moving inward 
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FIG. 4: Temperature at the neutrinosphere versus time for M = 1.55M0, (left upper figure), 
M = 1.7 Mq, (right upper figure), M = I.SMq, (left lower figure), M = l.QM©, (right lower 
figure) . 

because the matter is falling in. The star is shrinking until t ~ 0.16 ms, which is half of 
the oscillation period, and the density at the neutrinosphere is also at its minimum. Also at 
t ~ 0.16 ms, the outward heat pulse meets the infalling material density minimum. This can 
also be understood from the definition of the radius of the neutrinosphere, which is defined 
as the optical depth equal to unity. Therefore, when the temperature at the neutrinosphere 
is maximum, the corresponding density must be at the minimum value. 

While the resolutions we used in the simulations (which are limited by the computational 
resource) are good enough to model the global dynamics of the star accurately, we note that 
this is not the case near the stellar surface where the density is very small and changing 
rapidly. In particular, the grid resolution near the neutrinosphere is not very high. We have 
compared the numerical results obtained with a few different resolutions in order to examine 
the effects of resolution. 
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FIG. 5: Density at the neutrinosphere versus time for M = I.SSMq, (left upper figure), M = 
1.7 Mq, (right upper figure), M = 1.8Mq, (left lower figure), M = I.OMq, (right lower figure). 

In Figs. [71 we show the evolutions of the temperature and density at the neutrinosphere 
for the collapse model with M = 1.75Mq. The figures show that the different resolution 
results agree quite well qualitatively during about the first 1.5 ms. In particular, the period 
of the pulses does not depend strongly on the resolution. The maximum variation in period 
for different grid size is given by 0.28km/Cs, where Cg ~ lO^cm/s is sound speed near the 
surface, and it gives the maximum shift in period ~ 0.03ms. 

IV. EMISSION OF NEUTRINOS AND e± PAIRS 



A. Neutrino Luminosity 



The neutrino luminosity is given by 49], 
L,, = Attt'^c ■ 



27rr J 1 + exp(E^ - liy/kTy) ' 



(16) 
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FIG. 6: Temperature evolution (figure on the left) and density evolution (figure on the right). 
The arrow indicates the position at which the temperature (density) change occurs at the neutri- 
nosphere. 





FIG. 7: Evolutions of the temperature (figure on the left) and density (figure on the right) at the 
neutrinosphere for M = 1.75Mq, with three different grid resolutions. 

where /i^, is the neutrino chemical potential. Taking /ii, — 0, the neutrino luminosity emitted 
from the neutrinosphere is given by Lj, = iTrRlacT^/ld, where a — Aa/cis the radiation con- 
stant, (7 is the Stefan-Boltzmann constant, and 71, is the temperature of the neutrinosphere. 
If we assume equal luminosities for neutrinos and antineutrinos, the combined luminosity 
for a single neutrino flavor is 



L,j 77 



7 

Ltj = -nRlacT^. 
o 



(17) 



The effect of coherent forward scattering must be taken into account when considering 
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the oscillations of neutrinos traveling through matter (author?) [50|]. Although different 
flavor neutrinos have different R^y , yet they have approximately the same value of luminosity 

reviews and in depth discussions of t he p resent status of 



neutrino oscillations and their astrophysical implications see 



53 



54 



55 



, |56|] ) . Therefore the 



total luminosity is around three times of a single neutrino flavor luminosity 



L — Ly^^y^ -\- Ly^^y^ -\- Ly^^y^ 
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-TcRyacTj^ . 



(18) 



Using Ry and obtained from the last Section, we compute the neutrino luminosity as 
a function of time. The results are shown in Fig. [HI 
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FIG. 8: Neutrino luminosity versus time for M = I.SSMq, (left upper figure), M = 1.7Mq, (right 
upper figure), M = I.SMq, (left lower figure), M = l.OM©, (right lower figure). 



The peak luminosities range from 10^^ to 10^^ ergs/s; the pulsating period of the lumi- 
nosity is the same as that of the temperature and of the density. 
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B. Pair Production Rate 



Neutrinos and antineutrinos can become electron and positron pairs via the neutrino and 
antineutrino annihilation process (z/ + z/ — »• e~ + e"*"). The total neutrino and antineutrino 



annihilation rate can be given as follows 57], |58| 

= ^-^j^^D[kUt)f Q{r)4nr'dr (19) 

where 9(r) = 27r^(l — + 4x + 5)/3, x = ^1 — Rl/r'^, T^if) is the temperature at 

the neutrinosphere at time t, G\ = 5.29 x 10~^^ is the Fermi constant, ( is the Riemann 
zeta function, and D is a numerical value depending on the pair creation processes (e.g. 
experimental results indicate that Di = 1.23 for Uei'e and D2 = 0.814 for z/^ z/^ and z/^z/^). 
To obtain the total neutrino annihilation rate from all species, z/g, z/^ z/^ and Uf, we sum 
up the energy rate for each single flavor, 

= ''ffif\ D,+2D,)Rl(knr (21) 

Fig. [9] shows the rate of energy carried away by the electron/positron pairs produced 
through neutrino annihilation, which varies from ~ lO^^ergs/s to ~ lO^^ergs/s. 

It is interesting to note that almost all neutrinos can be annihilated into electron-positron 
pairs at the peak because of the extremely high density and high energy of the neutrinos. 
In particular the rest mass of the electrons/positrons is much smaller than kT^. 



V. MASS EJECTION AND ACCELERATION 



In order to calculate the mass ejected from the stellar surface by neutrinos/antineutrinos 
and pairs, a very detailed knowledge of the mass distribution near the surface is required. 
It was pointed out that the density profile in the crust plays a very important role in 
determining how much mass can be ejected [59!]. They use a static star model and an 
assumed simple power law density profile to demonstrate that the mass ejection can be 
significantly different. In our computer capability the minimum spacial grid size that can be 
achieved is 0.28 km. In order to estimate a precise location we choose to use the Piecewise 
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0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 

Time / ms Time / ms 



FIG. 9: Electron-positron luminosity versus time for M = 1.55Mq, (left upper figure), M = I.IMq, 
(right upper figure), M = 1.8Mq, (left lower figure), M = I.OM©, (right lower figure). 

Cubic Hermite Interpolating Polynomial (PCHIP) method to interpolate the density and 
the temperature data along the grids. The PCHIP method can provide a more accurate 
representation of the physical reality [goI. As compared with other interpolation methods 
(e.g., cubic spline data interpolation), the curve produced by the PCHIP method does not 
contain extraneous "bumps" or "wiggles", meaning that it could preserve the shape of the 
density and of the temperature profile, even when they change dramatically. However, it 
is unavoidable that even if we choose the best possible method, the true location might be 
slightly different from the real one. In Fig. [TD] we compare the original time evolution of 
the position of the neutrinosphere (-Rj,), temperature at Ry, density at Rj^, and the neutrino 
luminosity with the results obtained by using the PCHIP method. 

We can see that these quantities are qualitatively the same. We believe that the real 
continuous mass distributions near the stellar surface can be approximated by a continuous 
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FIG. 10: Results with PCHIP method and without data interpolation for M=1.55M0. 

mass distribution function, obtained from the numerical simulated data by using a PCHIP 
method. 

The mass ejection from a newly born quark stars was calculated in 59|]. They argue 



that neutrino-electron scattering is the most dominated process to deposit the neutrino 
energy in the crust. In this paper we argue that the dominated energy deposition process 
is the neutrino-antineutrino annihilation process. It should be noted that the optical depth 
(r), which is defined as r = nal, is the most important quantity to determine the energy 
deposition instead of cross section (cr) alone. Here n is the scattered particle density and / 
is the characteristic length. We can only eject mass above neutrinosphere, which is located 
only several grid sizes below the stellar surface when the neutrino luminosity is maximum. 
From Fig. 4 we can see that the density is several lO^^g/cm^ when the neutrino luminosity is 
maximum, where most of the mass is ejected. The density has dropped a factor of 10 in one 
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grid size from the neutrinosphere to the stellar surface. Therefore the scale length of density 
is about half of grid size, i.e. / ~ 0.14km. In a neutron rich matter, we can take electron 
fraction as 0.2 then we obtain nj ~ lO^^cm"^. On the other hand, the neutrino density is 
uniform from the neutrinosphere to the surface of the star and the density of neutrino is 
given by rij, = llaT'^/AkT ~ 10^^ {kT /15 MeV)^, where a is the Stefan-Boltzmann constant 
and k is the Boltzmann constant. At the neutrino luminosity maximum, kT ~15 MeV and 
the distance from the neutrinosphere to the stellar surface is several grid sizes, which is 
~lkm. Then n^l ~ lO^^cm^^. Since neutrino-antineutrino annihilation cross section (cf. 



Eq. 2 in [61]) and neutrino-electron scattering cross section (cf. Eq. 7 in 59|) are almost the 
same so we will ignore this process in calculating the mass ejection. 

A. Energy deposition in the crust and mass ejection 

Although most of the neutrinos/antineutrinos created above the neutrinosphere can es- 
cape, part of them can still be absorbed in the crust. We can estimate the amount of 
neutrino energy deposited in the crust due to the absorption in the following way. If we 
define Ru as i?NS > -Rm > Ru, then the absorbed neutrino energy onto the surface mass 
layer between Rm and i?NS could be expressed as 

E.{Rm) = / [l - e""^''"^ L{Ru) dt, (22) 

where t{Rm) = dr Kaffir) is the optical depth at i?M, L{Rm) = 2l7ri?|j a cT(i?M)'^/8 is 
the neutrino luminosity above -Rm and T{Rm) is the temperature at -Rm- Notice that we 
only performed data output for every 20 iterations in our simulations, corresponding to a 
time interval AT = 0.0075 ms. Hence, the time interval dt in the integral is taken to be 
dt = AT. 

The annihilated pairs created in the crust will be absorbed because of the much stronger 
interaction matter than that of neutrinos, and the pair energy {Ei±) deposited in the crust 
is given by 

Ei± = f QiRu) dt, (23) 



^(^m) = ^ ^,1^^ \ D, + 2D,){knf / Qir)Anr'dr. (24) 
Note that we only integrate r from -Rm to -Rns instead of integrating from -R,^ to oo. 
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The gravitational binding energy of the surface mass is 

GMAm , . 

= —B , (25) 

-Km 

where M = 4 vr p(r) dr, and Am{Ru) =4-7? /r^^ r"^ p{r) dr, respectively. 

Since the neutrino and pair absorption inside the neutron star actually happen simulta- 
neously, we combine the absorbed neutrino and pair energy together to be E'^^^"^^'^'^ 

^absorbed _ ^ _ ^26) 

As long as E^^^°^^'^'^{Rm) > Eg{Rm), the surface layer of Am^Ru) could be ejected from the 
neutron star. Hence from this criteria we can obtain the maximum ejected mass. 



B. Acceleration by pairs 

As we have mentioned in Section 4.2, the neutrino and antineutrino annihilation is very 
high at the peak of the neutrino pulses due to the extremely high density and the high 
energy of neutrinos. In fact most pairs are created outside the star, and therefore after the 
matter is ejected, it will be accelerated by absorbing the pairs created by the annihilation 
processes. The annihilation energy created from the neutron star surface Rns to r > Rns 
is given by Ei± = J Q^p^e±{r,t)dt, where 

^(^' = + 2^2) [kUt)]' f e(r') AT^r'^dr'. (27) 

Ic^n^ Jrns 

In the following we will briefly describe how the pairs accelerate the ejected matter. Before 
we present our calculations, we first describe the continuous mass ejection processes. The 
time slice interval of our output data is AT = 0.0075 ms. We can calculate the maximum 
amount of ejected mass only time slice by time slice. 

At Ti (Fig. [n] upper left), a layer of mass AM{Rm) is ejected when E''^'''"^'"^{Rm) > 
Eg{Rai)- The outer surface Rmij of the ejected mass is approximately Rns, and the velocity 
of the mass at the outer surface is almost c; the inner surface of the ejected mass is Rmis, 
where the velocity of the mass at the inner surface is the escaping velocity, which is almost 
half of the speed of light. 

T2 is the next time step when another layer of mass could be ejected. Before T2 (Fig. [11] 
upper right), when the mass layer ejected at Ti is flying outwards from the star at Ti+t < T2, 
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T = ^2 T = dT < Tj 



FIG. 11: Schematic illustration of mass ejection from the stellar surface. 

the distance to the ejecta can be approximated as RNsi+ct and t is the elapsing time counted 
from Ti. In the space between Rnsi and RNsi+ct, pairs are created and they can move faster 
than the ejected mass layer. Eventually they are absorbed by the ejecta and accelerates the 
mass layer in front. 

Another mass layer is ejected at T2 (Fig. [TT] lower left). When this ejecta moves outward, 
it will absorb pairs created between the stellar surface and this ejecta. However pairs created 
in the space between these two ejecta can still be absorbed by the first ejecta (cf. Fig. [TT] 
lower right). The total pair energy for accelerating the mass layer ejected at Ti is 

Too fTi+AT 

E'^''= E / QTArur2)dt, (28) 

Ti=Ti "^^^ 

where 

QrAn, r2) = —^^^iD, + 2D,){mf Q{r) Anr'dr, (29) 
is always calculated between the stellar surface and the ejecta or between two ejecta. 
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Pulse 


r 


Energy(ergs) 


Mass (g) 


Time (s) 


1 


17 


2.2 X 


1q46 


1.5 


X 


10^4 


0.16 


2 


6.4 X 10^ 


6.5 X 


IqA5 


1.1 


X 


1020 


0.50 


3 


1.6 


4.7 X 


1q48 


9.1 


X 


10^7 


0.83 


4 


1.6 


3.7 X 


1049 


7.1 


X 


1028 


1.18 


5 


1 X 10^ 


4.2 X 


1046 


4.6 


X 


10^9 


1.50 


5 


1.2 X 10^ 


2.7 X 


1048 


2.4 


X 


102^ 


1.51 


6 


2.8 X 10^ 


4.4 X 


1048 


1.7 


X 


1023 


1.85 


6 


1.4 


3.9 X 


1049 


1.1 


X 


1029 


1.86 


7 


3.8 X 10^ 


4.4 X 


10^1 


1.3 


X 


1028 


2.21 


7 


1.2 


2.4 X 


1049 


1.6 


X 


1029 


2.23 


8 


1.2 


8.9 X 


1048 


5.0 


X 


1028 


2.55 


9 


1.3 


3.7 X 


1048 


1.3 


X 


1028 


2.62 


10 


2.5 X 10^ 


2.3 X 


1048 


1.1 


X 


1022 


2.90 


10 


8.7 


2.6 X 


1050 


3.8 


X 


1028 


2.91 



TABLE I: Properties of the ejected mass. Ejected mass in some pulses arc divided 
into two part: first part with sufficient mass and move faster; the later part moves 
slower and cannot collide with first part. (M = 1.55Mo, rn=1.85 and 5^/^=160 
MeV). 

The Lorentz factor of each ejecta can be estimated as 

^pair _j_ ^absorbed 

r = + 1. 



(30) 



Wc summarize the results of the ejected masses in Tables 1-4. We can see that sometimes 
more than one layer of mass can be ejected in one pulse. 



VI. POSSIBLE CONNECTION WITH GAMMA-RAY BURSTS 
A. Duration of mass ejection 



We have pointed out that our numerical simulations are quite limited by the unphysical 
numerical damping. However, in realistic situations, the oscillations triggered by the collapse 
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Pulse 


r 


Energy(ergs) 


Mass (g) 


Time (s) 


1 


21 


1.6 


X 


10^5 


8.7 


X 


1022 


0.15 


2 


2.6 X 10^ 


1.0 


X 


1q45 


4.4 


X 


1021 


0.47 


3 


3.4 X 10^ 


3.2 


X 


10^5 


1.1 


X 


1022 


0.79 


4 


1.3 X 10^ 


5.3 


X 


10^5 


4.5 


X 


1022 


1.10 


5 


2 X lO'' 


3.6 


X 


1048 


2.0 


X 


1023 


1.41 


6 


2.2 


6.9 


X 


1049 


6.2 


X 


1028 


1.74 


7 


4.8 


2.0 


X 


1050 


6.0 


X 


1028 


1.89 


8 


2.3 X 10^ 


3.1 


X 


1048 


1.5 


X 


1025 


1.91 


9 


1.6 X 10*^ 


2.9 


X 


1048 


2.0 


X 


1021 


2.08 


9 


5 


4.7 


X 


1050 


1.3 


X 


1029 


2.09 


10 


1.1 


4.2 


X 


1048 


3.4 


X 


1028 


2.13 


11 


8 


1.2 


X 


1046 


1.9 


X 


1024 


2.41 


12 


31 


2.5 


X 


1044 


9.2 


X 


1021 


2.73 



TABLE II: Properties of the ejected mass. (M = 1.7M0, r„=1.85 and 51/^^=160 
MeV). 

must also be damped by some physical dissipation mechanisms. In this section, we would 
like to estimate how long the oscillations could last. 

First, we note that hydrodynamics effects (e.g., shock waves, mass-shedding etc.) which 



can be modeled by our simulations certainly would play a role in the damping |17l . |18 |. 
In particular, the study of [is!] suggests that the damping timescale due to hydrodynamics 
effects seen in their simulations is typically a few tens to hundreds ms. Another damping 
effect which exists in the case of rotational collapse is that due to gravitational radiation 
back-reaction. However, the damping timescale of this process is much longer than that of 
the hydrodynamics effects. The gravitational radiation damping is negligible and is in fact 
not taken into account in the works of jl?!. fisj]. There are still other damping mechanisms. 
It was first pointed out in \^2\ that the dissipation due to nonleptonic reaction is of great 
importance and the stellar pulsations of the quark stars would be strongly damped via 
s + u u + d in a few milliseconds. On the other hand, it was shown in (author?) 63 1 



and 



641 ] that in the high-temperature limit, which is exactly our case, the bulk viscosity is 
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Pulse 



1 

2 
3 
4 
5 
6 
7 
7 



9 

10 



6.1 X 10^ 
30 
7.8 

7.9 X lO'' 
5.9 
14 

2.2 X 10*^ 
8.2 

3.1 X 10^ 

3.9 
1.5 X 10^ 

59 



Energy(ergS; 



9.1 X 10^^ 

2.5 X 10^6 
1.3 X 10^*^ 

1.6 X 10^^ 

4.1 X 10^7 
3.5 X lO^s 
2.5 X 10^^ 

4.2 X 10^° 

4.2 X 10^*^ 
1.9 X 10^° 
1.0 X 10^^ 

7.3 X 10^3 



Mass (g) 



1.7 X 10^9 
9.5 X 10^3 

2.1 X 10^'' 

2.2 X lO^i 

9.3 X 10^5 
3.1 X 10^6 

1.3 X 10^1 

6.4 X 10^8 

1.5 X lO^i 
7.4 X 10^8 

7.3 X 10^1 

1.4 X 10^1 



Time (s) 



0.14 
0.46 
0.76 
1.06 
1.37 
1.68 
2.01 
2.01 
2.22 
2.22 
2.31 
2.64 



TABLE III: Properties of the ejected mass. (M = l.SM©, r„=1.85 and ^1/^=160 
MeV). 

dramatically reduced. If we take ~ 140 MeV, (p) ~ 10^^ g/cm^ and (T) ~ 50 MeV, the 
damping time scale is ~ 10 s [631]. 

If all of the above physical mechanisms cannot efficiently damp out the oscillations, then 
the pulse neutrino emission and the produced mass ejection considered in this paper would 
be the only processes to damp out the oscillations. In Tables 1-4 we can see that the energy 
carried away by the ejecta in first 3 ms is of the order of ~ 10^° ergs. The escaped neutrinos 
will carry away comparable amount of energy. The total oscillation energy is of the order of 
AEg ~ GM'^AR/ R^, where AR is the change of radius before and after the phase transition. 
For the models presented in this paper, AR/ R is of the order of 10%, which gives oscillation 
energy of the order of several 10^^ ergs. In other words the oscillations cannot last longer 
than several hundreds milliseconds even the neutrino emission and mass ejection are the 
only mechanisms to damp out the oscillations. 



29 



Pulse 


r 


Energy(ergs) 


Mass (g) 


Time (s) 


1 


2.6 X 10^ 


6.1 


X 


1q46 


2.6 


X 


1023 


0.13 


2 


6.4 


2.5 


X 


10^7 


5.2 


X 


1025 


0.43 


3 


6.3 X lO'' 


1.4 


X 


1046 


2.4 


X 


1020 


0.73 


4 


1 X 10^ 


2.4 


X 


10^7 


2.7 


X 


1024 


1.03 


5 


1.42 


1.9 


X 


1049 


4.9 


X 


1028 


1.32 


6 


4.3 X 10^ 


4.9 


X 


1048 


1.3 


X 


1025 


1.62 


6 


2.2 


7.0 


X 


1049 


6.5 


X 


1028 


1.63 


7 


1 X 10^ 


3.6 


X 


1049 


4.0 


X 


1023 


1.74 


8 


9.6 X lO'' 


5.2 


X 


1048 


6.0 


X 


1022 


1.94 


8 


76 


3.7 


X 


10^1 


5.6 


X 


1028 


1.95 


9 


8.8 X 10^ 


5.7 


X 


1048 


7.2 


X 


1023 


2.25 


9 


2.5 


1.9 


X 


1050 


1.4 


X 


1029 


2.25 


10 


7.5 


2.7 


X 


10^7 


4.6 


X 


1025 


2.55 


11 


1.2 


3.9 


X 


1048 


2.0 


X 


1028 


2.86 



TABLE IV: Properties of the ejected mass. (M = 1.9M0, r„=1.85 and B^/'^=160 
MeV). 

B. Short GRBs 



There are two kinds of GRBs, i.e., loner GRBs with duration larger than ~ 2 s and short 

n n 

GRBs with duration less than ~ 2 s [2, l3|. The isotropic 7-ray energy released by a short 
GRB is usually in the range of 10^^ — 10^^ ergs, i.e., about two or three orders of magnitude 
less than that of long GRBs 65] . Assuming that no more than ~ 10% of the kinetic energy 
can be converted to 7-ray radiation, then an amount of kinetic energy up to lO^'' — 10^^ ergs 



should be produced by t 
release is GRB 051221A 



le central engine. An example that has a relatively large energy 



661]. The isotropic 7-ray energy is 1.5 x 10^^ ergs, and an isotropic 



kinetic energy of 8.4 x 10^^ ergs has been estimated for this event 



66| 



It is widely believed that long GRBs may originate from the collapse of massive stars 



while short GRBs may be connected with the merger of binary compact stars 65|, 



67l |. However, it is still possible that some GRBs may be produced by other mechanisms. 
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For example, GRB 060614, a special nearby long GRB that is not associated with any 
supernovae, may be produced by an intermediate mass black hole that captured and tidally 



disrupted a star 



68| . Another interesting kind of central en gine mechanism involves the 



16 



69, 



70 



7l|. Since the phase 



phase transition of normal neutron stars to strange stars 
transition may be processed in a detonative mode, the details in this process are still largely 
uncertain and need further investigation. Especially, numerical simulations are necessary to 
help to understand the process. 



1. Basic equations 



The observations of X-ray, optical and radio after 



lave proved their cosmological origin 



iii,y 



23,1781,179, 



80 



81 



3, 



72 



73 



ows from some well-localized GRBs 



The so-called fireball model 



74, \L 



can basically explain the observational facts well, and thus it 
is strongly favored, and widely accepted today. In this model, the central engine gives birth 
to some energetic ejecta intermittently, like a geyser, producing a series of ultra-relativistic 
shells. The shells collide with each other at a radius of Ri^ and produce strong internal 
shocks. The highly variable 7-ray emission in the main burst phase of GRBs should be 
produced by these internal shocks. After the main burst phase, the shells merge into one 
main shell and continue to expand outward. It sweeps up circum-burst medium, being 
decelerated and producing external shocks. The observed long-lasting and steadily decaying 
afterglows (with a much smoother light curve as compared with the 7-ray light curve) should 
be due to these external shocks. 

After examining the numerical results of our simulations, we believe that the gravitational 
collapse of a neutron star induced by the phase transition from normal nuclear matter to 
quark matter can be an ideal mechanism for producing GRBs. In this Section, we give some 
detailed explanations. 

Let us first consider the general case of the collision between two typical shells. Following 
the description of [83|, we assume that the first shell is ejected with a mass of Mi, a bulk 
velocity of (3i and a bulk Lorentz factor of 71. The second shell is assumed to be ejected 
after a time interval of At, with the mass, velocity and Lorentz factor being M2, /32, and 72, 



respectively. Here, 71 = (1 — [3l\ 



^1/2 



72 



;i - and they satisfy 72 > 7i > 1. The 



distance between the two shells is then initially cAt. Since the second shell moves faster, it 
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will finally catch up with the first shell at a radius of 83 1 



,A,.^.i2l24,A, (31) 

p2 -Pi P2- Pi 72 - 7l 



This is the place where the internal shock appears and the GRB takes place. At this point, 
the total elapsed time measured in the static burster frame is 



. |. . ^ . A, (32) 

PiC p2 -Pi 72 - 7l 



However, due to relativistic effect, the observed elapsed time since the beginning of the phase 
transition is only tin/(27^) ^ 7|At/(7f — 7^). 

After the collision, the two shells will merge into one single shell and move at a new 
bulk velocity of /5bui (correspondingly, with the Lorentz factor 7bui)- During the process, a 
portion of the initial bulk kinetic energy will be dissipated as random internal energy. We 
denote the average Lorentz factor of the random internal energy as 71. Since momentum 



and energy are conserved in the collision, we have |83| . 



Mi7i/3i + M272/32 = (Ml + M2)7i7bui/3bui, (33) 

Mi7i + M272 = (Ml + M2)7i7bui. (34) 
It is then easy to get the solution for /3bui, 7bui and 7i as 

Mi7i/ji + M2J2P2 

Pbul = 77 — TT , (35) 

Mi7i + M272 



7bul 



Mi7i + M272 



Ml + Ml + 2MiM27i72(l - P1P2 



I M171 + M272 
Ml/71 + M2/72 



(36) 



M171 + M272 _ VMi-fi + M2I2 ■ \jMi/-ii + M2/72 
~ (Ml + M2)7bui M1 + M2 ■ ^ ^ 

The efficiency of transferring bulk kinetic energy into internal energy is 

e = (7i-l)/7i=l-7r'- (38) 



When two shells collide, the emission from the shock-accelerated electrons will correspond 
to a single pulse in the light curve of the GRB. The rising time of the pulse is mainly 
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determined by the time needed for the shocks to cross the two shells. Since the merged shell 
is moving toward us ultra- relativistically with a Lorentz factor of 7bui, only a small part of 
the shell (with an opening angle of ~ l/7bui) will be seen by us. The decay time of the 
pulse is then mainly determined by the arrival time lag of a photon emitted at the angle of 
~ l/7bui as compared with the photon emitted simultaneously at the line of sight. Usually, 



the decay time is longer than 



the width of the pulse, i.e. 



;he rising time, so we can use the decay time to characterize 



84 



85| 



'^puls 



^in/(27bulC). 



(39) 



2. GRBs resulting from phase transition? 

According to our simulations, at least more than 10 shells can be ejected during the phase 
transition of a neutron star in the first 3 ms. Among these shells, a few will have isotropic 
energies larger than 10^^ ergs and expand ultra- relativistically with Lorentz factors larger 
than > 100 (cf. Tables 1-4). Note that although our simulations only last for about 3 ms, 
the actual duration of the oscillation process may last for several hundreds milliseconds, as 
argued at the end of our Section 6.1 . During this period, tens or even hundreds of ultra- 
relativistic shells might be ejected, each with an energy larger than ~ 10^^ ergs. The total 
energy enclosed in these relativistic shells may be 10^'^ — 10^^ ergs. We suggest that the 
collision between these shells can give birth to a GRB. 

To describe this process in a quantitative way, we first study the collision between two 
specific (but typical) shells. We assume Mi = M2,7i = 300,72 = 271 = 600. We further 
assume that the second shell is ejected after a time interval of At = 1 ms. It is then 
straightforward to find that they will collide at a radius of i?in ~ 8'jfcAt/3 ~ 7.2 x 10^^ cm, 
and at the time of tin ~ 8jfAt/3 ~ 240 s. Note that due to relativistic effect, the observed 
elapsed time since the beginning of the phase transition is only tin/(27^) ~ 4At/3 ~ 1.3 ms. 
The merged shell will move at a Lorentz factor of 7bui ~ V^li ~ 420, with the co-moving 
internal energy characterized by 7i ^ 3/2^/2 ^ 1.06. The efficiency of transferring the bulk 
kinetic energy into internal energy is e = 1 — 1/1.06 ~ 5.6%. The pulse will have a width of 
Tpuise ~ 2At/3 ^ 0.6 ms. 

In realistic cases, the shells are ejected with variable masses, variable Lorentz factors, 
and variable time intervals, respectively. The conditions then become very complicated. 
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For example, if 72 is much larger than 71 (72 S> 71, but we still assume M2 = Mi), then 
the solution can be expressed as Rin ~ 2j'lcAt, tin ~ 27^ At, 7bui = ^^7172, 7i = ^72/471, 
Tpuise = At7i/72. In this case, the pulse width will be very small. On the other hand, if 72 
is only slightly larger than 71 (still with M2 = Mi), then 7bui ~ 7i,7i ~ 1,-Rin and tin will 
be very large so that a pulse much wider than At can be generated. However, in this case, 
the energy transfer efficiency will be extremely low (e ^ 1) for this single pulse. 

In the standard fireball model, it is generally believed that the total duration of a GRB 
refiects the active period of the central engine. In our simulations, the shells are ejected 
mainly in a few hundred milliseconds. So, this mechanism is most proper for explaining 
short GRBs. The duration may reasonably range between a few ms and several hundred 
ms if the neutrino emission and the mass ejection are the only processes to damp out the 
stellar oscillations. It is also possible for the duration to extent to more than 1 s. For 
example, if we take 71 = 350, 72 = 400 and At = 300 ms, then the observed elapsed time for 
this pulse will be larger than tin/(27^) = 7|/(7|-7i)At ^ 4.3At ~ 1.3 s. Note that the total 
isotropic energy of the shells can be lO^'' — 10^^ ergs, and this is also consistent with the 
energy requirement of most short GRBs. If some extent of anisotropy exists in the process 
(which is quite possible if the effects of strong magnetic field of the compact star are further 
included, see the discussion in the last paragraph of this section), then the energy release 
can even meet the requirement of those rigorous energetic events, such as GRB 051221A 
Q^. Also, as noted above, if two shells are ejected with similar Lorentz factors, they may 
collide at a very late time. We further notice from Tables 1 — 4 that there are many high 
energy shells that moves at trans-relativistic speeds (with Lorentz factors significantly less 
than ~ 10). Late collisions can also be produced by these energetic shells when they finally 
catch up with the decelerating external shock (that produces GRB afterglow). Such late 
collisions may manifest themselves as fiaring activities (emerged in the afterglow phase), as 
frequently observed 1000 — 10000 s after the trigger of GRBs [86|, l87 |. 

In our simulations, the shells are ejected periodically, with a period of ~ 0.3 ms. However, 
it is quite unlikely that we could observe any obvious periodicity in the 7-ray light curve. 
The reason is that the shells have different masses, velocities, and energies. The periodicity 
will then be most likely smeared out at the time of collision. 

The energy releases in our simulations are basically isotropic. The reason is that we did 
not include the effects of rotation and macroscopical magnetic field in our calculations. When 
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these effects are considered, the energy releases should show some features of anisotropy. 
This is an interesting point that needs to be investigated in further details. 



VII. DISCUSSIONS 

In this paper we have studied the possible consequences of the phase-induced collapse of 
neutron stars to hybrid stars. We have found that both the density and the temperature 
inside the star will oscillate with the same period, but almost 180° out of phase, which 
will result in the emission of intense pulsating neutrinos. The temperatures at the peaks 
of pulsating neutrino fluxes are 10-20 MeV, which are 2-3 times higher than non-oscillating 
case. Since the electron/positron pair creation rate sensitively depends on the temperature 
at neutrinosphere, the efficiency of pair creation increases dramatically. We want to point 
out that the intense pulse neutrino and pair luminosity can be maintained due to the os- 
cillatory fluid motion, which can carry thermal energy directly from the stellar core to the 
surface. This process can replenish the energy loss of neutrino emission much quicker than 
the neutrino diffusion process. Part of neutrino energy, roughly (1-1/e), and pairs inside 
the star will be absorbed by the matter very near the stellar surface. When this amount 
of energy exceeds the gravitational binding energy of matter, some mass near the stellar 
surface will be ejected, and this mass will be further accelerated by absorbing pairs created 
from the neutrino and antineutrino annihilation processes outside the star. Unlike inside the 
star, the surface properties, e.g. position of the neutrinosphere, surface temperature etc., 
are very sensitive to the surface perturbation. The neutrino and pair luminosities can be 
varied from pulse to pulse. This results in the large variation of the Lorentz factor of the 
mass ejecta. The internal collisions among these mass ejecta may produce the short time 
variabilities of the Gamma-ray Bursts, which can be as short as submilliseconds. Although 
mass ejecta are ejected periodically, each ejecta can have different masses and Lorentz fac- 
tors as explained before. Therefore, the intrinsic period could not be observed. Although 
we can only simulate the oscillating stars for a few millisecond, we can speculate that this 
may be a possible mechanism for short Gamma-ray Bursts based on the following reasons. 



891] have estimated that the viscous damping time for such oscillating system is of order of 
10 s. By assuming that if neutrinos and pair emissions are the only damping mechanisms, 
the pulsations can last less than ~ 3s [89], which is roughly the characteristic time scale of 
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short Gamma-ray Bursts. 



, with the 



The phase-transition from a neutron star to a strange star was simulated in 
conclusion that this process is most likely not a gamma-ray burst mechanism. They mimic 
the phase-transition by the arbitrary motion of a piston deep within the star, and they 
have found that the mechanic wave will eject ~ IO^^Mq baryons, which causes the baryon 
contamination for the gamma-ray bursts. In our simulations, we assume a sudden change of 
equation of state to mimic the phase-transition, and we use the Newtonian hydrodynamic 
code to study the response of the stellar interior, after such a sudden change of the EOS. In 
our simulations we find that the mass ejection by the motion of the fluid is very small. We 
estimate that the major mass ejection would result from the heating of neutrinos and pairs 
on the stellar crust, which is not modeled in the simulations. Our total energy output and 

(cf. Tables 1-4). However, the neutrino energy 



total mass in ejecta are close to that of 



injection is pulsating, and hence the mass ejection is also pulsating. The masses of individual 
ejecta range from ~ 10~^Mq to ~ 1O~^M0, with output energy in the range of 10^^ ergs 
to 10^° ergs. Therefore, some ejecta cannot be relativistic, and they cannot contribute to 
GRBs. However, there are still many relativistic ejecta in each simulation model, which can 
have Lorentz factors >100, and with a total energy of ~ 10^° — 10^^ ergs. We conclude 
that this could be a possible mechanism for short GRBs. 

Finally, we want to remark that our numerical simulations describe a spherically symmet- 
ric, non-rotating, and collapsing stellar object. Also the effect of the magnetic fields was not 
taken into account. Therefore, the radiation emission produced in this model is isotropic. 
However, a realistic neutron star should have finite angular momentum and strong magnetic 
field, and hence these two factors could produce asymmetric mass ejection. This effect will 
be considered in future work. 
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